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The electromagnetic form factors of the deuteron, particularly the quadrupole form factor, are 
studied with a help of a phenomenological Lagrangian approach where the vertex of the deuteron- 
proton-neutron with D-state contribution is explicitly taken into account. The result shows the 
importance of this contribution to the quadrupole form factor in the approach. 
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I. INTRODUCTION 

The study of electromagnetic form factors of nucleon and light nuclei, like deuteron and He-3, are crucial for the 
understanding of the nucleon structures. Deuteron, as the most simplest nuclei, has been a subject of many years 
(for some recent reviews [H-Q). Since it is a weekly bound state of the proton and neutron, the study of the deuteron 
can shed light on the study of the nucleon as well as on the nuclear effects. Moreover, as a spin-1 particle, the 
deuteron structures are different from the spin-1/2 nucleon and He-3, and from the spinless pion meson. There are 
many discussions on the deuteron structures, like its wave functions, binding energy, the electromagnetic form factors, 
and the parton distributions, in the literature. Those works are usually based on the phenomenological potential 
models with quark, meson, and nucleon degrees of freedom and based on some effective field theories etc. [H-ll^. The 
realistic deuteron wave functions, with the help of meson exchange potential model, have been explicitly given by 
Refs. [Il-[3. 

In our previous works [il[i3: a phenomenological Lagrangian approach is applied for the electromagnetic form 
factors of the deuteron, where it is regarded as a loosely bound state of a proton and a neutron, and the two constituents 
are in relative S-wave for simplicity. The coupling of the deuteron to its two composite particles is determined by the 
known compositeness condition from Weinberg , Salam [I^ and others [1^ [2l| . Our phenomenological effective 
Lagrangian approach has been proven to be successful in the study the weekly bound state problems, like the new 
resonances of Y(3872), and A+(2940), the EM form factors of pion as well as some other observables 1^. 

It should be stressed that since only one-body S— wave operator contribution is considered in our previous study 
the estimated quadrupole moment of the deuteron is negligibly small when compared to data. According the 
non-relativistic potential model calculation [l^, one sees that the deuteron quadrupole moment is very sensitive to 
the D— wave component of the deuteron. Therefore, S— state contribution is not sufficient. In order to avoid the 
discrepancy, several two-body arbitrary and phenomenological Lagrangians were introduced, by hand, to compensate 
the discrepancy (T^ . 

The purpose of this work is to re-study the deuteron electromagnetic form form factors with the phenomenological 
approach. Here both the S— and D— state contributions to the vertex of the deuteron-proton-neutron are simultane¬ 
ously taken into. It is expected that by considering the D— state contribution in the vertex, the estimated deuteron 
quadrupole could be sizeably improved. This paper is organized as follows. Section 2 briefly shows our theoretical 
framework, particularly the D— state contribution to the vertex. Numerical results and some discussions are given 
in section 3. 


II. FRAMEWORK OF THE APPROACH 

Deuteron, as a spin-1 particle, has three independent form factors. The matrix element for electron-deuteron (ED) 
elastic scattering, as shown in Fig. 1, can be written as 


( 1 ) 
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under the one-photon exchange approximation. In eq. (1) k and k' are the four-momenta of initial and final electrons 
and J7^(P, P') stand for the deuteron EM current. Its general form is 

J^{P,P')= (2) 

- (^Gi me'* ■ e - ^^e ■ qe'* ■ (P -h P'), 

-G2(Q')(^e^e'*-9-e';e-g^ , 


where is the deuteron mass, e(e') and P{P') are polarization and four-momentum of the initial (final) deuteron, 
and = —q^ is momentum transfer square with q = P' — P. The three EM form factors Gi_ 2,3 of the deuteron are 
related to the charge Gc, magnetic Gm, and quadrupole Gq form factors by 

Gc = Gi + , Gm = G 2 , (3) 

Gq = Gi-G2 + (l + r)G3, 

with the factor of t = Q^/AM^. These three form factors are normalized at zero recoil {Q^ = 0) as 

Gc(0) = 1, Gq{ 0) = MjQd = 25M, (4) 

Gm(0) = ^fid = 1.7U, 


where M^v is the nucleon mass, Qd and fid are the quadrupole and magnetic moments of the deuteron. 

The unpolarized differential cross section for the eD elastic scattering can be expressed by the two structure 
functions, A{Q‘^) and as 


da 

c?n 


= O'M 


+ P(Q^)tan^ 


( 5 ) 


where ctm = cos^(0/2)/[4P^ sin^(0/2)] is the Mott cross section for point-like particle, E and E' are the incident 

and final electron energies, 9 is the electron scattering angle, = 4PP'sin^(0/2), and a = e^/Air — 1/137 is the 
fine-structure constant. The two form factors A{Q^) and B{Q^) are related to the three EM form factors of the 
deuteron as 

Am = Gim + Ir^Gim + IrGl^m 

Bm = \r{l + T)Gl,m- ( 6 ) 

Clearly, the three form factors Gc,m,q cannot be simply determined by measuring the unpolarized elastic eD differ¬ 
ential cross section. To uniquely determine the three form factors of the deuteron one additional polarization variable 
is necessary. For example, one may take the polarization of T 20 [l| ■ 

Take an assumption that the deuteron as a hadronic molecule-a weakly bound state of the proton and neutron, 
one may simply write a phenomenological effective Lagrangian of the deuteron and its two constituents of the proton 
and neutron as 


Cd{x) = gnOm J dy^Diy^)pix + y/2)T'^Gn^{x - y/2) + H.c., 


( 7 ) 
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where is the deuteron field, Cn^{x) = n^^ix), and C = 17 ^ 7 ° denotes the matrix of charge conjugation, and x is 
the centre-of-mass (C. M.) coordinate. In the above equation is the vertex for the deuteron-proton-neutron and 
the correlation function characterizes the finite size of the deuteron as a pn bound state. The correlation 

function $_d(j/^) depends on the relative Jacobi coordinate y. 

If only the S— wave contribution is considered, the simplest form of the vertex is F^ ^ 7 ^ which has been employed 
before [2. When both the S— and D— states contributions are considered, then the vertex becomes more complicated. 
According to the work of Blankenbecler, Gloderber, and Halpern the vertex of the deuteron-proton-neutron is 

F^ = F^'^+F^'^ (8) 


where the first and second terms stands for the contribution from S— and D— states, respectively. They are 



(9) 


and 


p2,/x _ 
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( 10 ) 


with p being a measure of the U—state admixture, k is the relative momentum between the proton and neutron, and 
= M;^d with 6 being the binding energy of the deuteron as shown in Fig. 2. Here, it should be mentioned that in 
the rest frame of the deuteron, the non-relativistic reduction gives 


T-il ,2 /-y ^ -* 


ea2 



( 11 ) 


It means a combination of two spin-1/2 states, proton and neutron, forms a spin triplet state. Similarly, in the 
non-relativistic limit, ^ 7 ^ — means the proton and neutron couple to spin triplet state and this spin triplet 

state re-couples {k) to form a state with the same quantum numbers of the deuteron. 


The coupling of the deuteron to its two constitutes, in eq. (7) is determined by the known compositeness 
condition Z — Q proposed by Weinberg, Salam and others [ISl - l^ . This condition implies that the probability to 
find a proton and neutron system inside the deuteron is unity. Thus, the coupling oi go is determined according to 
Zd = 1 — = 0 , with 


^'0(^0 ) — 9d^'d±{^d) ( 12 ) 

being the derivative of the transverse part of the mass operator (see Fig. 2). Usually, the mass operator splits into the 
transverse and longitudinal parts ofI]p^(fc) = g“^S£)_L(fc^)-|-(fc^), with —and 5 “^A:„ = 0 . 

We see that the coupling of the deuteron to its constituents of the proton and neutron, 5 ^), is well determined by the 
compositeness condition. 

A basic requirement for the choice of an explicit form of this correlation function is that its Fourier transform vanishes 
sufficiently fast in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. Usually 
a Gaussian-type function is selected as the correlation for simplicity. One may choose d>D(fc^) = exp{— A:|/A 2 ) for 
the Fourier transform of the correlation function, where kE is the Euclidean Jacobi relative momentum and A stands 
for the free size-parameter which represents the distribution of the constituents in the deuteron. 
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Here the analytical expression for the coupling is 

(13) 

9d 

where ^ and 2 stand for the derivatives of the transverse parts of the mass operator from the contributions 
of the S— and D— states, respectively. The explicit expressions are 
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where pLN,d = ^/A| and 
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with 


i ' ' i 

B'{a, jd) = ^ {a T /3) + Aa(d 


Z^ = 


IT aso ^ ^ ^ 

—- 2 - 1 - 0 ! + ^, 


(17) 


and asD = A|/A|j. Here we simply ignore the p^- dependent term since p is expected to be small, and we consider 
that the correlation functions of the S— and D— states are not necessarily the same, therefore we have totally 
three parameters Ag, Ad and p in this calculation. Then, we can calculate the matrix element of photon-deuteron 
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FIG. 4: Estimated deuteron charge form factor Gc(Q^). The solid and dotted curves are the results of our calculations and of 
the phenomenological parametrization . The data are open circle f^. open square[3^, open diamond [U], plus [ 3 ^. triangle 
upd^, filled circle [ 3 ^ . and filled squared^, respectively. 


interaction as shown in Fig. 3. We have 




E E 

{N=p,n) (ij'=l,2) 
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where the photon-nucleon current of 


1 

= PiMQ")r + 


(18) 


(19) 


is employed with and F 2 ,Ar being the known nucleon Dirac and Pauli form factors and N = p,n, stand for the 
proton and neutron, respectively. 


III. NUMERICAL RESULTS AND DISCUSSIONS 

We calculate the matrix element of eq. (18) and consider the one-photon exchange approximation for the photon- 
deuteron current as shown eq. (1). Thus we can get the corresponding deuteron three form factors Gi_ 2,3 as well as the 
deuteron charge GdQl, magnetic Gm{QI and quadrupole Gq(Q^) form factors. There are so me p arameterizations 
for the nucleon form factors of F’i, 2 (Q^) in the literature for the proton and neutron by |25l - l^ . In the present 
calculation, we employ the ones of Blunden The three mo del-dependent parameters, A 5 = 0.10 GeU, Ap = 
0.08 GeU and p = 0.03, are fixed by fitting to the experimental data. The obtained charge, magnetic and quadrupole 
form factors are shown in Figs. (4-6). 


It should be stressed that, in this work according to the discussions of Ref. [l^, we explicitly include the D— 
state contribution to the deuteron-proton-neutron vertex as shown in eq. (10). Comparing to our previous work 
in [l^ . we found that this contribution is very important for the understanding of the quadrupole moment and 
quadrupole form factors. The estimated Gm( 0) and Gq( 0) are about 1.53 and 21.38, respectively. These two values 
are reasonable comparing to the normalization conditions of 1.714, and 25.83 given in eq. (4). If we only take the 
S— wave contribution into account, we hardly reproduce the experimental measurement for the quadrupole moment 
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FIG. 5: Estimated deuteron magnetic form factor Gm{Q^)- The solid and dotted curves are the results of our calculations and 
of the phenomenological parametrizationf^ . The data are circle [ 3 ^ , sauaref37|. diamond [3^. and triangle [3^. respectively. 


at the zero-recoil limit, although the estimated magnetic moment is consistent with the data. Here, the negligibly 
small value of the quadrupole moment in Ref. [l^ is improved due to the inclusion of the D— state contribution. 
Meanwhile, the charge and magnetic moments also remain reasonably. 

In summary, we explicitly consider, in this work, the D— state contribution to the vertex of the deuteron-proton- 
neutron, as well as the S— wave one simultaneously, and find that our four-dimensional phenomenological Lagrangian 
approach can reasonably reproduce the deuteron charge, magnetic, particularly, quadrupole form factors simulta¬ 
neously. The estimated quadrupole moment is much improved due to the inclusion of the D— state contribution. 
It should be stressed that our present approach is a fully relativistic and it is different from the potential model 
calculations based on the three-dimensional framework. 

Of course, the present calculation can be further improved, since we still cannot reproduce correctly the crossing 
point of the charge and magnetic form factors of the deuteron as discussed in Ref. [ 2 ^ . It is found that the experimental 
data for Gc and Gm show the existence of a zero, for = 0.7 GeV^ and — 2 GeV^, respectively. This is 
probably due to the fact that our selected correlation functions are still simple. Moreover, the explicit form of the In¬ 
state contribution, as shown in eq. (10), is not unique M- A more sophisticated calculation is in progress. Finally, 
it is expected that the future calculation of the deuteron generalized parton distribution functions with help of this 
approach is promising. 
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FIG. 6: Estimated deuteron quadmpole form factor Gq{Q^). Notations are the same as Fig.4. 
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